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^ ■ Abstract 



We study the lattice gauge model proposed recently by Kazakov and Migdal for 
inducing QCD. We discuss an extra local Zn which is a symmetry of the model and 
propose of how to construct observables. We discuss the role of the large- N phase 
transition which should occur before the one associated with the continuum limit 
in order that the model describes continuum QCD. We formulate the mean field 
approach to study the large- iV phase transition for an arbitrary potential and show 
that no first order phase transition occurs for the quadratic potential. 
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1 Introduction 



Recently Kazakov and Migdal have proposed a very interesting lattice gauge model 
for inducing QCD. The model is defined by the partition function 

Z = /n^(*)nd*(*)e^^ (LI) 

X,[l X 

where the field <E>(x) takes values in the adjoint representation of the gauge group SU(N) 
and the link variable U^{x) belongs to the group. The potential V[$] can be expanded 
at small lattice spacing a in the power series 

V[$\ =m^$ 2 + A$ 4 + ... (1.2) 

with m being the bare mass of the field <£>. A strong coupling solution to this model 
has been constructed by Migdal at iV = oo and then extended to next orders of the 
1/iV-expansion ||. 

To obtain the solution two assumptions have been made. The first one is that the 
critical point = D (D is the dimensionality of euclidean space-time), which separates 
the strong coupling region, is associated with continuum QCD. The second one is that 
the path integral over $(x) is saturated by a single x- independent configuration — the 
master field. 

On our mind along solving the Kazakov-Migdal model, one should better understand 
qualitative properties of the model like how to construct observables or how to recover the 
area law in the continuum. As usual the question of which phase transitions occur in the 
lattice system is of special importance. The knowledge in advance of the phase structure 
might help to understand what kind of the distribution of eigenvalues of the matrices 
$ij{x) would correspond to a physical solution. In this paper we apply a standard lattice 
gauge theory technique to study phase transitions in the Kazakov-Migdal model. 



A subtle point with the model ( PL1D is that it possesses an extra local Zn symmetry 



U^x) -> Z^x)U^x). (1.3) 

The important role of this symmetry in the Kazakov-Migdal model has been pointed 
out recently by Kogan, Semenoff and Weiss [Q. These authors argued that due to the 
^Tv-symmetry the model undergoes a phase transition which was assumed to occur si- 
multaneously with the transition from strong to weak coupling solution. This phase 
transition separates two phases which differ by the nature of confinement — the so-called 
local confinement versus usual area law. 

The main subject of the present paper is a study of properties of the Kazakov-Migdal 
model which come from the presence of the Z/v-symmetry. To calculate observables, in 
particularly the string tension, we propose a general procedure based on the Wilson loops 
in the adjoint representation SU(N) which are invariant under the Zjy-symmetry and 
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can be expressed entirely in terms of averages over $(x). We show that, being proper 
normalized to be of order 1 in the naive continuum limit, these quantities are ~ 1/N 2 
in the strong coupling expansion. For the quadratic potential V[<&], we look for the 
corresponding critical point at iV = oo using the mean field method. We show that no 
phase transition occurs in D = 4 up to the point = D where the model becomes 
unstable. 

We discuss that, if the large- iV phase transition occurs before the one associated with 
the continuum limit, a consistent picture of inducing continuum QCD by the Kazakov- 
Migdal model emerges. In this case one gets below the large- N phase transition a theory 
with usual confinement (area law) rather than local confinement as in each order of the 
strong coupling expansion. We formulate the mean field approach to study the large- N 
phase transition for the case of an arbitrary V[<&]. 

We speculate that the phase transition associated with the continuum limit should be 
identified with the one separating confinement and Higgs phases and discuss the general 
scenario of inducing QCD by the Kazakov-Migdal model. We point out that the mean 
field analyses of the Higgs phase transition is reduced to a one-matrix problem which is 
solvable for cubic and logarithmic potentials. 

2 Zn and observables 

The local Zn symmetry^ has far-reaching consequences for the model ( |1 . 1| ) . First of 
all, the average of any non-invariant quantity like the Wilson loop in the fundamental 
representation vanishes except for a loop passing the same contour back and forth (i.e. 
with vanishing minimal area). Moreover, this property holds independently of how many 
phase transitions the system undergoes under the way to the continuum since the local 
symmetry can not be broken spontaneously. Therefore, the average of the fundamental 
Wilson loop always vanishes. 

To obtain the continuum Wilson loops from the Kazakov-Migdal model, we propose 
to consider averages of the adjoint Wilson loops 

W A (C) = ^(\trU(C)\ 2 -l)) (2.1) 

which are invariant under the Z^. We use the normalization that W A ~ 1 in the naive 
continuum limit. 

For the sake of simplicity let us postpone for a moment discussing properties of the 
adjoint Wilson loop for the model (|1.1|) and review some results || for the case of the 

2 Since the transformation can be done independently on each link of the lattice, such a symmetry was 
called in Ref.g that of the 3-rd kind. 
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pure adjoint single-plaquette lattice action 



^4 = -^ I 2 (2.2) 



□ 



where Pa ~ 1 as N — > oo in order to have a nontrivial large- iV limit. Due to the 
factorization at large N, the following formula holds 

W A {C) = (^trU(C) ) 2 + 0(N- 2 ) (2.3) 

where the average of the Wilson loop in the fundamental representation on the r.h.s. 
should be calculated for the Wilson action 

S Wilson = -NpJ2^ u ( n ) (2-4) 

□ 

with the coupling (3 given by the self-consistency equation 

P = /3 A W{n;P) (2.5) 

where W(n-/3) is the plaquette average for the Wilson action (with the coupling being 
P). When Eq. (|2.5|) possesses a nontrivial solution which is valid in the weak coupling 
region so that asymptotically 

P — ► Pa~^ as P A -> oo, (2.6) 

one sees from Eq. (|2.3| ) that the adjoint Wilson loop display the area law at N = oo with 
the adjoint string tension being 

Ka = 2K Fundamental- (2-7) 

The perimeter law which is expected for the adjoint Wilson loop at finite N enters the 
term of order 0(N~ 2 ). 

Let us now return to the model (|1.1|). Our idea is to define at iV = oo the 'fundamental' 
Wilson loop, which enters, say, the correlator of electromagnetic currents represented via 
sum over paths or determines the string tension, by Eq. ( |2.3| ) taking the square root 
of the adjoint Wilson loop. This procedure is unambiguous since the imaginary part of 
the fundamental Wilson loop never shows up to the leading order of 1/iV-expansion. It 
is crucial for this procedure the large- N phase transition to occur before the continuum 
limit sets in. Only in this case the induced continuum theory would possess normal area 
law while otherwise one gets local confinement. 

Thus, we can define the set of QCD observables at the kinematical level despite the 
Z/v-symmetry while this procedure works only at iV = oo. It is a dynamical question 
whether the large- N phase transition makes this construction sensible. 
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3 Large- N phase transition 



Let us start again from the single-plaquette action (|2.2|) . Since the only solution of 
Eq. ( |2.5| ) in the strong coupling expansion is (3 — 0, the adjoint Wilson loop Wa{C) 
vanishes in this region according to Eq. ( |2.3| ). One can directly verify this result estimating 
the order in 1/N of the strong coupling expansion which gives Wa(C) ~ 1/N 2 for the 
action (|2~2l ). Therefore, as we concluded in Ref. f|, there should be a first order phase 



transition for the action (|2.2| ) at some value f3 A < 2. 

This phase transition was observed by the Monte-Carlo simulations for iV = 2 4- 6 
H with the critical value of f3* A growing from 1.6 for AT = 3 to 1.8 for iV = 6 what 
is close to the estimate (3* A < 2. More detailed studies of Eq. (|2.5| ) were performed in 
Ref. 0. An uncertainty of the value of (5* A is related, in particular, to an ambiguity of 
criterion for the first order phase transition at large N. The standard criterion stating 
that the phase transition occurs when free energies of two phases coincide (the 'Maxwell 
rule') may not work at A" = oo because the free energy itself is ~ A^ 2 and a barrier 
which separates the two phases becomes infinite in this case. An alternative criterion says 
that the phase transition is associated with the point where a metastable weak coupling 
solution terminates. Using the 'Maxwell rule', Samuel obtained (3 A = 1.54 0. Analogous 
studies of this phase transition using mean field / variational technique yields (3* A = 2.8 
while the criterion based on terminating the metastable region gives (5 A = 1.7 in better 
agreement with the Monte-Carlo data. The large- A" phase transition is not related to 
breaking of a symmetry and was shown to be associated with dynamics of Z^-monopoles 
which are condensed in the strong coupling phase ||. 

An occurrence of the analogous phase transition for the Kazakov-Migdal model has 
been advocated by Kogan, Semenoff and Weiss ||. We support the consideration of 
this paper by the following argument. Let us estimate the order in 1/N of the adjoint 
Wilson loops ( |2.1| ) in the strong coupling expansion of the model ( Jl.lj) with the quadratic 
potential V[<E>] using the induced action [l|] 

SUU] = -Ij: 1 -^4^. (3.1) 

It is easy to see now using the standard lattice technique that W A {C) ~ 1/A^ 2 in the strong 
coupling expansion while it is of order 1 in the naive continuum limit. This estimate is 
quite similar to the one discussed in the previous section for the case of the single-plaquette 
adjoint action. Thus, we conclude the model undergoes a first order phase transition with 
decreasing toq. 

Let us now estimate the location of this phase transition. Following the scenario of 
Ref. 0, we start from the induced action (|3~T| ) and decrease mj. A very interesting 
question is whether the large-A" phase transition occurs before the point tuq = D after 
which the quadratic action in Eq. ( |1 . 1|) is not bounded from below. The point is that the 
action (|3.1|) is not, say, the classical action for an external field problem. One should inte- 



5 



grate exp(—Si n d[U]) over U^{x) with the Haar measure according to the definition (|L1 



As nig — > oo, typical configurations of U^x) are uniformly distributed on the group 
being independent on each link of the lattice. This disordering suppresses the contribution 
of long loops to the sum over paths on the r.h.s. of Eq. ( |3.1| ) — the longer the loop the 
stronger the suppression. In usual lattice gauge theory the L/-matrices become ordered 
under the way to the continuum limit. The ordering is, however, only for distances smaller 
than the correlation length while for large distances the disorder is needed for confinement. 
This ordering occurs either gradually or is enhanced by the presence of a first-order phase 
transition. It is, however, a dynamical question whether this phase transition occurs in 
the region ttIq > D. 

A simplest logical possibility for the model (|1.1|) might be that the large-iV phase 
transition occurs at some value of m\ > D which is large enough in order that the action 
( p.l| ) can be approximated by the single-plaquette term, i.e. coincides with the action 
(P3D with 

^ = sib- (3 ' 2) 

This would mean that the Kazakov-Migdal model were induce a lattice gauge theory with 
the action ( j2.2|) . Substituting the above numerical value [3* A ~ 2, one would get ml ~ 0.6 
which is too small. The next terms on the r.h.s. of Eq. (|3.1| ) are essential for such m\ so 
that this situation is excluded by dynamics. 



4 The mean field analyses 



To study the large- N phase transition, we apply to our problem the mean field method, 
which usually works pretty well for first order phase transitions, quite analogously it was 
applied to the single-plaquette adjoint action in Ref. [Q. 

To construct the mean field, we use the variational approach which was advocated in 
the context of modern field theory by Sakita [ECO]. Let us introduce the trial partition 
function 

"\{dU^x)e b ^^ tvU ^ x)? (4.1) 

which is a product of one-link integrals. We have chosen for them a simplest form which 
possesses the symmetry. The Jensen's inequality yields then the following bound on 



the partition function fll.l ) 



Z>Z r[[d$(x)e ^ x 21 rtJl i/o, (4.2) 

J X 

where (. . .)o means averaging w.r.t. the trial action. Since the exponent contains the sum 
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of one-link averages, it can be expressed via the one-matrix integral 



2 JdUe^ trU \ ^\tiU 



V 2 = — » ' ' , (4.3) 

JdUe^ trU ^ 

by the formula 

(^tr^(x)U IJi (x)^(x + fj)U^(x) ) = r/ 2 tr $(jc)$(ac + fi) . (4.4) 

The idea of the variational mean field method is to fix 6 a from the condition that the 
trial ansatz (|4. 1|) would give the best approximation to Z in the given class. Calculat- 
ing the derivative w.r.t. 6 a and taking into account that rj depends on bA according to 
Eq. (fOp, one finds the maximum of the r.h.s. of Eq. (|4.2|) at 



6 A = 7 v • (4.5) 

This is the final formula that relates 6a to the potential V[$] providing the r.h.s. of 
Eq. ( [4.3| ) is known as a function of 6a- Its standard interpretation is that one obtains the 
one matrix integral on the r.h.s. of Eq. (|4.3| ) replacing the matrix [U^x)]^ in the partition 
function ( |1 . 1| ) by the mean field value rjdij on each link of the lattice except the given one 
(0, + fi) while Eq. ( |4.3| ) gives a self-consistency condition at this link. 

The one-matrix integral on the r.h.s. of Eq. Q4.3|) was first calculated by Chen, Tan and 



Zheng H (see also Ref. J7|). In the weak coupling region where &A>2orl/2<?y<l, 
the result can be represented as 

2( v -r] 2 ) = ^. (4.6) 

OA 

We have written the self-consistency condition in the form which would be convenient to 
study the phase transition. 

Let us consider the case of the quadratic potential V[$] when the mean field analyses 
is drastically simplified. The gaussian integral on the r.h.s. of Eq. (|4.5|) can easily be 
calculated to give 

b A = —,\ dae'^ I°- l (a) L» (4.7) 

T] z JO 

where 1(a) are modified Bessel function. 



Eqs. (|4.3| ), ([LTD can be analyzed similarly to those of Ref. ||. The r.h.s. of Eq. 
multiplied by rj 2 monotonically increases with decreasing rn^/rj 2 with maximal value ~ 0.1 
at ml/r] 2 = 4 and then diverges. This value is too small to provide a solution to Eq. (|4.6|) 
with rj > 1/2 as it should be for the weak coupling solution. 
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We conclude, therefore, that for the quadratic potential V[<&] there is no first order 
large-N phase transition for m\ > 4 when the model is well-defined. Notice that we ex- 
clude a possibility that the large- A r phase transition occurs exactly at ml = 4. According 
to Eq. ( |4.(j| ) this would occur if were > 1/2 at this point which is not the case. 

Let us mention that the gaussian integral over $ on the r.h.s. of Eq. ( |4.5| ) for the 
quadratic V[$] can be represented in the 1/wjq expansion as the sum over paths: 




where the sum goes over the open loops with the endpoints x and x + /j,. This sum is not 
exactly the same as what would appear if the mean field method were applied directly to 
the action ( |3.1| ). The point is that if a link is passed back and forth, it does not contribute 
to the r.h.s. of Eq. ( |3.1|) because U is unitary while it is taken into account on the r.h.s. 
of Eq. ( |4.8D . In particularly first few terms of the series expansions of r.h.s. of Eq. (|4.7|) 
in D = 4 read 




which differs to order 0{l/mfy from what one would obtain from Eq. ( |3.2| ). All the extra 
terms are, however, positive so that fe^ can be considered as an upper limit which is 
enough for our purposes. 



5 Discussion 

The main conclusion from the fact that the large- N phase transition does not occur for 
the model ( |1 . 1|) with the quadratic potential V[<&] is that it remains in the strong coupling 
phase with local confinement and can not induce QCD. This phase is pretty trivial: the 
averages of all the Wilson loop vanish at N = 00 except of those with vanishing minimal 
area, in contrast to the strong coupling expansion of the lattice gauge theory with the 
Wislon actionP]. Therefore, one should incorporate the self-interaction of $ and look 
for the large-iV phase transition. The mean field method which leads in the case of an 
arbitrary potential V[<&] to Eqs. ( |4.5|) and ( |4.6|) can be used to estimate the location of 
the large- iV phase transition. 

As we discussed already, the large-A r phase transition should occur before the one 
associated with the continuum limit in order that the Kazakov-Migdal model induce 
normally confining QCD. One might tempt to relate the latter phase transition with the 

3 The existence of a master field for the Wilson action in the strong coupling region was first advocated 



by Kazakov, Kozhamkulov and Migdal 1 1 1 



one separating confining and Higgs phases which always occurs in gauge theories. A 
continuum theory with confinement can be usually obtained by approaching this phase 
transition from above while approaching from below one gets the deconfining Higgs phase. 
If one accepts this scenario, the main problem with inducing QCD by the Kazakov-Migdal 
model is to study whether the large- iV phase transition is indeed separated from the Higgs 
one. If the two phase transitions coincide, this would mean that one passes from the 
phase with local confinement directly to the Higgs phase while the phase with normal 
confinement were missing. Monte-Carlo simulations of the Kazakov-Migdal model might 
help to answer this question. 

The mean field method could be useful to study the Higgs phase transition as well. 
In the Higgs phase one would get a condensate of the $- field, with a nonsymmetric 
distribution of eigenvalues of $* which violates the SU(N) quite similarly, say, to the 
Higgs phase transition in the Georgi-Glashow model. One should add then one more 
self-consistency condition to determine $*. Substituting $(x) by an (independent on x) 
mean-field value <3>* for all sites of the lattice except given one, we get the self-consistency 
condition 

J_ tr 1 = J E C 5 i) 

N X-Q^ JdQeNtrl-Vffl+DrP&S,*) ' ^ > 

This condition means that the saddle-point configuration of the integral over $ coincides 
with and A plays the role of a spectral parameter. 

The matrix integral in Eq. ( |5.1| ) coincides with the partition function of the hermitian 
one-matrix model in an external field, While its solution for a quartic potential V[$], 
which has been taken into account in Ref . J2| , is not yet known, this model was explicitly 
solved in genus zero (the N = oo limit) for a cubic potential |12| and for a logarithmic 
potential [O]. We hope the obtained results could be useful in this context. 
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Added note 

When this paper was being prepared for publication, there appeared more papers |L4[ 
on the Kazakov-Migdal model. The results by Gross agree with ours for the quadratic 
potential while the Monte-Carlo study by Gocksch and Shen seems to indicate that for 
N = 2 the 'large- N 1 phase transition coincides with the Higgs one. 
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